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ABSTRACT 

Nearly all globular clusters (GCs) studied to date show evidence for multiple stellar populations, in stark 
contrast to the conventional view that GCs are a mono-metallic, coeval population of stars. This generic feature 
must therefore emerge naturally within massive star cluster formation. Building on earlier work, we propose 
a simple physical model for the early evolution (several 10 8 yr) of GCs. We consider the effects of stellar 
mass-loss, type II and prompt type la supernovae, ram pressure, and accretion from the ambient interstellar 
medium (ISM) on the development of a young GCs own gas reservoir. In our model, type II SNe from a 
first generation of star formation clears the GC of its initial gas reservoir. Over the next several 10 8 yr, mass 
lost from AGB stars and matter accreted from the ambient ISM collect at the center of the GC. This material 
must remain quite cool (T ~ 10 2 K), but does not catastrophically cool on a crossing time because of the high 
Lyman- Werner flux density in young GCs. The collection of gas within the GC must compete with ram pressure 
from the ambient ISM. After several 10 8 yr, the Lyman- Werner photon flux density drops by more than three 
orders of magnitude, allowing molecular hydrogen and then stars to form. After this second generation of 
star formation, type II SNe from the second generation and then prompt type la SNe associated with the first 
generation maintain a gas-free GC, thereby ending the cycle of star formation events. Our model makes clear 
predictions for the presence or absence of multiple stellar populations within GCs as a function of GC mass 
and formation environment. While providing a natural explanation for the approximately equal number of first 
and second generation stars in GCs, substantial accretion from the ambient ISM may produce fewer chemically 
peculiar second generation stars than are observed. Analyzing intermediate-age LMC clusters, we find for the 
first time evidence for a mass threshold of ~ 10 4 M Q below which LMC clusters appear to be truly coeval. This 
threshold mass is consistent with our predictions for the mass at which ram pressure is capable of clearing gas 
from clusters in the LMC at the present epoch. Recently, claims have been made that multiple populations 
within GCs require that GCs form at the center of their own dark matter halos. We argue that such a scenario is 
implausible. Observations of the young and intermediate-age clusters in the LMC and M31 will provide strong 
constraints on our proposed scenario. 
Subject headings: Galaxy: globular clusters — globular clusters: general — stars: evolution 



1 . AN OBSERVATIONAL PUZZLE 

Globular clusters (GCs) have historically been considered 
coeval, mono-metallic, gravitationally bound collections of 
stars. In the past several years, high precision photometric 
and spectroscopic observations have led to a radical revision 
of this picture. 

High resolution spectra of stars within nearly all GCs stud- 
ied to date reveal internal spreads in light element abundances 
such as C, N, Na, O, Mg, and Al, beyond what can be ex- 
plaine d by measurement errors (see review in iGratton et alj 
2004). The magnitude of the internal spread varies con- 
siderably from cluster to cluster, though there ar e notice- 
able trends with cluster ma ss and orbital properties dCarrettal 
l2006t ICarretta et ail 1201 Oct) . Intriguingly, internal spread in 
the Fe-peak elements and the type II supernovae product 
Ca is limited to o nly the most massive GCs (tt>Cen, M22 , 
M54, and Terzan 5;lMarino et al||2 009: Carr etta et al]|2009at 
Ferraro et al. 2009; Carretta et al. 2010a, c). The most massive 
clusters also show indirect evi dence for a large spread in He 
abundances (Gratton et al. 2010). Abun dance variations h ave 
been detected in main sequence stars (IGratton et al.ll200lb . in- 
dicating that the observed variation arises from stars forming 
out of different material, as opposed to being due to some 
unknown mixing process, which could only occur along the 
giant branch. While much attention has been paid to internal 
abundance spreads recently, such internal spreads have been 



known for over 30 years dCohenlll978tlKraftlll979l:lPetersonl 
[19801 IFreeman & Norrisll981l:ISmith & Norrisll98 2b. 1983). 

Photometry from the Hubble Space Telescope (HST) has 
demonstrated that many GCs contain multiple sub-giant 
branches and at least two (uiCen and NGC 2808) contain mul- 
tiple main sequences (see review in Piotto 2009). Variations 
in the light element abundances have been conclusively as- 
sociated with the multiple sequences ob served in the color- 
magnitude diagram (CMD; Yong et al. 2008; Marin oet alJ 
120081: ICarretta et ai]|2009 b; Milone et aL Wiq) . demonstrat- 
ing that these two phenomena are intimately linked. The pop- 
ulation with enhanced abundance patterns is more centrally 
conce ntrated than the 'nor mal' populat ion in the GCs NGC 
1851 (Zoccali et all [2009b and o;Cen ([S ollima eTail [20071) . 
and perhaps many others as well (Carretta et al. 2009b|K 

In the old Milky Way (MW) GCs the relative numbers of 
stars with 'anomalous' and normal abundances ratios is ap- 
proximately equal, with little dependence on m etallicity (e.g., 
ISmifh & Nornslll982bl; ICarretta et al.ll2009bl) . This fact im- 
poses strong constraints on formation scenarios, as we explain 
in later sections. The relative numbers of anomalous and nor- 
mal stars in GCs in other galaxies (e.g., the LMC, M31) is not 
known, but would provide new insight. 

Additional insight has come from the strong observed cor- 
relations between various light element abundances within in- 
dividual GCs. The most striking is the anti-correlation be- 
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tween Na and O (e .g., iKraft et al.1 fl993t llvans et alj QUI 
iCarretta eT al. 2009b). This correlation arises naturally when 
material of standard (e.g., solar) abundance ratios is mixed 
with material that has been processedj at temperatures > 
10 7 K. At such temperatures the Na-Ne and CNO cycles are 
active, with the former producing Na and the latter depleting 
O. At similar temperatures the Mg-Al cycle is active, explain- 
ing observed correlations between these elements as well. 

The discovery of significant Li and F abundances in the 
second stellar generation (i .e., asso ciated with stars of high 
Na and low O abundances; Pasquini et al. 2005; Smith et al. 
2005) has provided yet another puzzle and clue. These el- 
ements are fragile, especially Li which burns at T > 10 6 K. 
The existence of such fragile elements in the atmospheres 
of second generation stars that also show strong O deple- 
tion and Na enhancement suggests that these stars formed 
out of at least two kinds of material; i.e., from matter ex- 
posed to T > 10 7 K, and additional material that was never 
heated above T ~ 10 6 K. The interpretation of Li is however 
complicated by the possibility that Li may, under special cir- 
cumstances, actually be p roduced within AGB stars (see e.g., 
IVentura & D'Antonall2010h . 

Multiple stellar populations have also been detected 
in intermedia t e-age (~ 1 Gyr) and old LMC clusters 

f lackey et all [2 008; Goudfr ooii et alj 120091 iMilone et alj 
09; Mucciarelli et al. 2009), a nd possibly in the old GCs 
within the Fornax dSph galaxy (Letarte et al. 2006). These 
observations indicate that the multiple stellar population phe- 
nomenon is not specific to the MW. Observations of LMC 
clusters are of course hampered by the much larger dis- 
tance modulus to the LMC. Despite this limitation, the 
intermediate-age clusters offer a new window into the inter- 
nal age spreads because the main sequence turn-off point at 
these ages is a strong function of time. Small age differences 
are therefore readily noticeable in the CMD. Analysis of HST- 
based CMDs have shown that the spread observed in the main 
sequence turn-off of intermediate-age LMC clusters can be 
explained with an internal age spread of a few 10 8 yr. 

The only class of star clusters known not to contain mul- 
tiple stellar populations are the open clu sters in the MW 
dde Silva et alJl200lllMartell & Smithll2009l) . which have typ- 
ical masses of a few thousand solar masses, although some, 
such as NGC 7789, have masses of ~ 1O 4 M . 

These observations have led to the unavoidable conclusion 
that the majority of GCs studied to date harbor multiple stel- 
lar populations. For all but the most massive ones, GCs are 
still considered chemically homogeneous in Fe-peak and el- 
ements arising primarily from SNe type II such as Ca. The 
light element variations have been detected in both young and 
old cl usters, both metal-poor and metal-rich (Martell & Smith 
2009), and are noticeably absent in the open clusters. 

From these observations the following timeline in the early 
evolution of GCs has emerged. Within GCs, a first genera- 
tion of stars form. Type II SNe then remove any remaining 
gas from the GC. After the epoch of type II SNe, mass from 
evolved stars is cycled through temperatures of T > 10 K 
and then is returned to the gaseous reservoir of the GC. Af- 
ter a few 10 8 yr, a second generation of stars forms from a 
mix of processed and un-processed material. Star formation 

1 Throughout the text we will use 'processed' and 'un-processed' to refer 
specifically to the nuclear processing of material at the temperatures required 
to explain the observed light element abundance variations, i.e., at T > 10 7 K. 



permanently ceases after the formation of the second genera- 
tion. This process does not occur in open clusters, nor in the 
field. This basic timeline has, in one form or ano ther, been 
discussed by many authors (e.g., iCottrell & DaC osta 1981; 
[Smith|[T987tlCarretta et alJl20T0M) . 

Any theory of GC formation must be embedded into our 
broader cosmological theories of structure formation. Many 
lines of evidence suggest that galaxy formation is an hierar- 
chical process where dark matter halos serve as sites for as- 
sembling baryons and converting baryons into stars. One of 
the inevitable predictions of this bottom-up scenario is that 
MW GCs likely formed in environments very different from 
the ~ 200 kms -1 dark matter dominated halo w here th ey now 
reside. One possibility, which we discuss in J J2.1.11 is that 
GCs form in the centers of small dark matter halos. This sce- 
nario would imply that isolated GC are embedded in extended 
dark matter halos. This paper emphasizes another possibility: 
GCs form within small gas-rich dwarf galaxies — the build- 
ing blocks of the present MW. We emphasize that the MW is 
an evolving galactic system so that estimates of ram pressure 
and gaseous accretion must consider the likely environment 
in which a GC formed at z ~ 2- 10, rather than its current 
location in the MW today. 

We have outlined above only the most basic sketch of what 
must occur to explain the observations. In the next section 
we critically assess previous, more detailed scenarios for the 
development of multiple stellar populations within GCs. Fol- 
lowing this assessment, we describe our own model for the 
early evolution of GC stellar populations that includes several 
novel ingredients, and is, at least in certain respects, more 
plausible than other scenarios. We also present an analysis of 
intermediate-age clusters in the LMC that provides confirma- 
tion of a key aspect of our model. We conclude by comment- 
ing on various observations that may shed new light on this 
exciting observational puzzle. 

2. PREVIOUS EFFORTS 

There are at least three major issues that require explanation 
before any satisfactory theory for the formation of multiple 
populations in GCs can be accepted. These are: 1) under- 
standing how the gaseous reservoir within a GC can remain 
within the shallow potential well for several 10 8 yr; 2) iden- 
tifying which stars are responsible for processing material at 
T > 10 7 K; and 3) explaining how a second generation of stars 
can form with a current total mass comparable to the first gen- 
eration. In this section we will assess the plausibility of pre- 
vious efforts at addressing these outstanding issues. 

2.1. Can mass lost from first generation stars remain bound 
to the GC? 

A serious challenge to any theory for the formation of mul- 
tiple stellar populations is the shallow potential wells of GCs. 
If gas is to remain bound within GCs, then one of the follow- 
ing two scenarios must occur: 1) the gas must both be ejected 
from stars below the GC escape speed and must remain cool 
so that the internal dispersion is < 10 kms" 1 ; or 2) the GC 
must be at the gravitational center of a more massive system 
with the escape speed from the larger system sufficiently large 
to retain stellar winds of a range of speeds. 

An additional obstacle is ram pressure. Many GCs cross the 
disk of the MW at high velocity, which results in a strong ram 
pressure force felt by gas within GCs. Under certain condi- 
tions, ram pressure may therefore prohibit the formation of a 
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gaseous reservoir w ithin GCs (see e.g., Frank & Gisler| [T976t 
iGnedin et alj|2002h . 

The shallow potential wells of GCs and their susceptibil- 
ity to ram pressure has led some to resurrect the notion that 
GCs form at the center of their own dark matter halos at 
high redshift (e.g., [B ekki & Norris 2006; Bekki et al.1 l2007t 
ICarretta eT al. 2010b). The feasibility of this scenario is as- 
sessed in the following section. 

2.1.1. Formation of GCs at the centers of small dark matter halos? 

Early theories regarding the formation of GCs included the 
possibility that they form within extended dark matter halos 
at high redshift (Peebles 1984). This scenario fell out of favor 
following the observation of thi n tidal tails surroun ding many 
GCs fe.g.. lGrillmair et al.lll995l:IOdenkirchen et al.ll2003l) . be- 
cause numerical simulations showed that su ch tidal tails d o 
not form if GCs reside within extended halos (lMoorefl l996). 

The formation of GCs at the center of their own dark mat- 
ter halos has received renewed interest thanks to numerical 
work that has shown that extended halos of GCs could be 
fidally stripped away by the pre sent epoch (iBromm & Clarkd 
120021: iMashchenko & SiTTsI 120051) . In this picture, the con- 
straint from tidal features demonstrates only that GCs at the 
present epoch are not embedded within extended massive ha- 
los, but offers no insight into the formation environment of 
GCs. 

This scenario has been adopted in order to explain the 
formation of m u ltiple stellar populations within GC s (e.g. , 
Freemanl [1991 iBekki & Nornst 120061 iBekki et alj 120071: 



Boker 2008; Carretta et al. 2010b) for the principle reason 



that GCs embedded within massive halos do not experience 
ram pressure and may also easily retain stellar winds. In 
our opinion, the motivation for dark matter halos from ram 
pressure arguments is in error. Recent arguments regard- 
ing ram pressure have been based directly on ram pressure 
calculations for GCs orbiting the MW at the present epoch. 
Indeed, no GCs today show evidence of a gaseous reser- 
voir, and ram pressure is thought to be resp onsible (see e.g., 
ITavler & Woodl[T975l IPrank & Gislerlll976l) . These calcula- 
tions assume mass-loss rates appropriate for ~ 10 10 yr old 
populations, and orbital velocities characteristic of the present 
MW. In contrast, at the epoch relevant for the formation of 
multiple populations (e.g., the formation epoch of GCs), the 
characteristic orbital velocities were almost certainly substan- 
tially smaller, probably by an order of magnitude, and the 
mass-loss rates were an order of magnitude larger. We will 
demonstrate quantitatively in a later section that these two 
effects combine to significantly reduce the efficiency of ram 
pressure so that all MW and LMC clusters that show evidence 
for multiple stellar populations were likely impervious to this 
effect in their formation environment, without appeal to dark 
matter halos. In short, ram pressure is not an important phys- 
ical process for the mass ranges and formation environments 
of the current set of clusters with known multiple populations. 

There are additional arguments against the formation of 
GCs within dark matter halos. The most striking is the obser- 
vation of multiple populations within intermediate-age LMC 
clusters. These clusters surely are not embedded within mas- 
sive halos, and yet they clearly are capable of forming mul- 
tiple populations. If these clusters can form multiple popu- 
lations without being embedded in a massive halo, why not 
MW GCs? 

If GCs formed at the center of their own halos, then the ex- 
istence of GCs within dwarf spheroidals would be extremely 



difficult to understand because they would have very short 
dynamical friction times. For example, Fornax contains five 
GCs. In the absence of an extended massive halo, the friction 
timescales are already uncomfortably short if For nax resides 
within a standard dark matter halo (several Gyr; Goerdt et al. 
2006). If these GCs were embedded in massive halos, their 
dynamical friction timescales would be <C 1 Gyr since fric- 
tion scales with the total GC mass. 

A final argument against GC formation within extended 
halos is the structure and kinematics of very isolated GCs 
both in the MW and M31. NGC 2419 is a remote MW GC 
at a distance of k, 90 kpc from the Galactic center, has a 
mass of ~ 10 6 M Q , and half-mass radius of rs 20 pc. Based 
on radial velocity measurements of stars within NGC 2419, 
Baumgardt et alT (120091) measured a mass-to-light ratio that 
was consistent with a pure stellar population. It is hard to 
imagine how tidal stripping could effect this cluster, and so the 
lack of any evidence for a massive dark matter halo for NGC 
2419 argues against GCs generically forming within massive 
halos. Recently, an extremely isolated GC has been found 
within M31, at a distanc e of rj 200 kpc from the center of 
M31 (Mackey et al. 2010). This object, far removed from any 
strong tidal fields, should be another ideal candidate to search 
for evidence of dark matter within GCs. 

The formation of GCs within massive halos does how- 
ever provide one significant advantage over models that do 
not invoke formation within halos. At high redshift (e.g., 
z= 10), the gas accretion rate onto a lOkms -1 halo from 
the intergalactic medium is of the order PigmV47t/? 2 w 6 x 
10~ 3 M Q yr" 1 . The accretion rates are much higher at z = 10 
than the present epoch because the density of the IGM is 
higher by a factor of (1 +z) 3 = 10 3 . After several 10 8 yr, this 
would result in a substantial accumulation of gas, which could 
help explain th e rel ative numbers of first and second genera- 
tion stars (see i j2.31 l and the high abundances of Li and F. 

In this section we have argued that 1) GCs need not form at 
the centers of their own halos in order to survive ram pressure 
stripping and 2) regardless of ram pressure considerations, 
not all GCs could have formed within their own halos given 
both dynamical constraints and observations of intermediate- 
age GCs in the LMC and M3 1 . Nonetheless, we emphasize 
the possibility that some GCs formed at the centers of their 
own halos. Nuclear star clusters reside at the centers of galax- 
ies, are massive (~ 10 6 - 10 8 M Q ), compact (half-light radii 
of a few pc), and span a wide range of ages (10 7 - 10 10 yr; 
iBoker et alj 2004; Walcher et all 12005. 200$. The youngest 
nuclear star clusters must certainly have formed where they 
are now observed. The formation environment of the older 
clusters is less clear because they may have migrated toward 
the center via dynamical friction. In any event, the existence 
of nuclear star clusters at least suggests that GC-like objects 
can in some cases form at the centers of their own dark matter 
halos. As we discuss below, the most massive GCs in the MW 
may be examples of such objects. 

2.2. Identifying the stellar type responsible for processing 
material to T > 10 7 A" 

Based on nucleosynthetic constraints, there are two plau- 
sible candidate donors of processed materia l: massive (> 
20M Q ) rotating stars dDecressin etalj 12007), and massive 
(w4- 8Afo) AGB stars ( Vent ura et alJl200U ID' Antona et alj 



120021) . Unfortunately, both of these proposed sites suffer from 
serious theoretical uncertainties (e.g., compare the AGB mod 
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els of iKarakas & Lattanziol d2007l) and IVentura & D'Antonal 
(2008tJ)! The relevant cross-sections for the production 
of Na are uncertain by a factor o f ~ 1000 (see discus- 
sion in IVentura & D'Anto na 2008b), making detailed com- 
parisons between models and data difficult. For example, 
Ventura & D'Antona (2008a) were able to construct AGB 
models whose ejecta abundance patterns were consistent with 
the abundances of observed second generation stars only after 
they modified several key yet uncertain nuclear cross sections. 
Models for massive rotating stars suffer from significant un- 
certainties not only in the mass loss rates but also in the effi- 
ciency of meridional circulations to mix processed material to 
the envelope. 

The massive rotating stars model suffers from several short- 
comings. First, the timescale for these stars to lose mass is 
comparable to the timescale for type II SNe. It is difficult to 
imagine how new stars could form in an environment where 
massive stars are constantly exploding. Even if this could oc- 
cur, the second generation would then show significant over 
abundances in a-elements compar ed to the firs t generation, 
but this is not observed (see also lRenzinill2008l who make a 
similar point). 

Moreover, the short timescales involved would imply that a 
different physical mechanism is at work in the intermediate- 
age LMC clusters, where CMD studies have shown internal 
age spreads of order 10 8 yr dMilone et al.ll2.Q09b . In contrast, 
the AGB stars evolve on timescales similar to the observed 
age spread in LMC clusters, providing a natural explanation 
for the observed spread. 

Finally, the short evolutionary timescales for the massive 
rotating stars makes it difficult to explain the observed Li and 
F abundances. The most plausible explanation for the rela- 
tively high abundances of these elements is that the second 
generation formed from a mix of un-processed and processed 
material. As we discuss in later sections, it is very difficult to 
accrete a sufficient amount of un-processed material in only 
~10 7 yr. _ 

Recently, de Min k et aT] ((2009) proposed a scenario 
whereby massive binaries shed a large fraction of their mass 
due to non-conservative mass and angular momentum trans- 
fer. We believe that this scenario suffers from many of the 
shortcomings of the winds from massive stars scenario, as 
outlined above, and is therefore also disfavored. 

We instead favor AGB stars as the source of processed ma- 
terial out of which the second generation stars form. Besides 
the fact that they are the only remaining option, the evolu- 
tionary timescale of massive AGB stars is comparable to the 
inferred internal age spread within intermediate age GCs in 
the LMC. 

2.3. Constraints imposed by the relative numbers of first and 
second generation stars 

As emphasized by many authors, it is very difficult to un- 
derstand the observational fact that within most GCs the num- 
ber of first and second generation stars is approximately equal. 
For example, assuming a standard Kroupa (2001) IMF, and 
assuming that AGB stars shed all of their mass aside from a 
0.6M© white dwarf remnant mass, only w 9% of the total ini- 
tial stellar mass of the system can be returned to the surround- 
ing GC from AGB winds with initial masses of 4- 8M(tQ. If 
100% of the AGB ejecta were converted into stars, the popu- 

2 Our returned mass fraction is substantially smaller than the fractions 
quoted in Bekki & Norris (2006) because they assume AGB stars down to 



lation ratio should be approximately ten first generation stars 
for every one second generation star. 

Several scenarios have been proposed to account for the 
discrepancy between the above expectation and the observed 
ratio. Early proposals focused on ad hoc variations to the 
IMF in order to produc e the observed relative numbers (e.g., 
ISmith & Norrisl Il982al) . Recent work has even attempted 
to constrain the IMF of the first generat ion in this way 
dD'Antona & C aloi 2004; Prantzos & Charbonnel 2006). 

More recently, iD'Ercole et al.l ((2008) proposed a solution 
that required the first generation to have been substantially 
more massive (by a factor of 10- 100) than its present mass. 
If the first generation was born with mass segregation already 
in place, then stellar evolutionary and dynamical effects could 
result in the loss (i.e., unbinding) of a significant number of 
first generation stars. The second generation, being born with 
higher concentration, would be relatively immune to such a 
process, and so the majority of second generation stars would 
remain bound. The preferential loss of first generation stars 
could then result in comparable numbers of first and second 
generation stars at late times. 

The proposal that second generation stars form purely from 
AGB ejecta suffers from several shortcomings, including the 
required level of fine-tuning. If we assume that the first gen- 
eration donates 10% of its mass to the GC gas reservoir via 
AGB winds and a star formation efficiency of 30%, the re- 
sulting total mass in the first and second generations will be 
equal if the first generation losses 97% of its initial mass. This 
assumes that no stars from the second generation are lost; any 
mass lost from the second generation would require an even 
larger fraction of ma ss-loss from the first generation. In the 
ID'Ercole et al.l d2008l) scenario, the first generation of every 
GC must lose > 97% of its mass so that the observed popula- 
tion ratios are always of order unity. Without a natural mecha- 
nism to explain why the first generation always sheds enough 
mass so that its final mass is comparable to the second gen- 
eration, the mass lost must be fine tuned ad hoc. One would 
also expect isolated GCs such as NGC 2419 that experience 
weak tidal fields to contain many more (by factors of 10- 100) 
first generation compared to second generation stars. Abun- 
dance measurements of stars within NGC 2419 currently do 
not exist but would clearly shed light on this issue. 

A second concern with the D'Ercole et al. solution is the 
large remnant population left behind by the disproportionately 
large fraction of massive first generation stars at the center of 
the GC. In this scenario, a GC with present mass 10 5 M Q and 
equal numbers of first and second generation stars would re- 
quire an initial stellar mass of ~ 10 7 M Q first generation stars 
in order to generate enough gas from AGB winds to produce 
a second generation with total mass ~ 10 5 M Q . For a standard 
Kroupa (2001) IMF, approximately 10% of the total initial 
mass is constituted by M > 8M Q stars. If 10% of these mas- 
sive first generation stars are concentrated toward the center 
because of mass segregation and then explode, a population 
of neutron stars and black holes will be left behind, and re- 
main bound, with a total mass of ~ 1O 4 M ! 

IMq can contribute to the GCs gaseous reservoir. The main sequence life- 
time of such stars is substantially longer than ~ 10 8 yr, and in addition it is 
believed that hot bottom burning, which is the source of many of the abun- 
dance anomalies, does not occur in such low mass stars. We therefore believe 
our fraction, which is based on the expected AGB mass range for hot bottom 
burning, to be more accurate. Expanding the adopted mass range to 3 — 8Mq 
results in a mass return fraction only slightly larger (13%) than what we quote 
for the 4 - 8Mq range. 
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Such a large population of massive remnants would not 
only be detectable as an additional dark mass component con- 
centrated toward the dynamical center, but would also imply 
a very high incidence of low-mass X-ray binaries (LMXB). 
Neither of these predictions are supported by the observa- 
tions (e.g., van der Marel 2004; Verbunt & Hut 1987), al- 
though there are notable caveats. The observed frequency of 
LMXBs can be naturally explained by the high rate of binary 
formation in GCs compar ed to the field, although this mech- 
anism is quite uncertain (Fabian et al. 1975; Verbunt & Hut 
1987). We have assumed that a significant fraction of neu- 
tron stars remain bound to the cluster, which may not be true 
if the majority of GC neutron stars are born with high kick 
velocities (as appears to be the case for pulsars in the MW; 
lHansen & PhinnevTT9 97). Appealing to IMF variations to ex- 
plain the relative numbers of first and second generation stars 
would also produce significantly more neutron stars and black 
holes, as pointed out by D'Antona & Ca loi (2004). A more 
careful analysis of the expected remnant population in GCs 
may yield st ronger co nstraints on thes e models. 

Recently, ISills & Glebbeekl (120101) have investigated the 
importance of stellar collisions in boosting the number of fast 
rotating massive stars and intermediate mass stars (i.e., AGB 
stars) in GCs. These authors find that stellar collisions are not 
sufficiently numerous for their products to help explain the 
relative number of first and second generation stars. 

In the following section we propose an alternative scenario 
for resolving this discrepancy that is based on adding a signif- 
icant amount of gas via accretion from the ambient ISM. 

3. A PLAUSIBLE MODEL FOR THE FORMATION OF MULTIPLE 
STELLAR POPULATIONS IN GLOBULAR CLUSTERS 

In the previous section we evaluated a number of current 
scenarios that address aspects of the formation of multiple 
stellar populations within GCs. In this section we propose 
an alternative scenario for the formation of multiple stellar 
populations in GCs. 

3.1. Model overview 

Our scenario can be summarized as follows. A first gen- 
eration of stars forms out of gas that has already been pre- 
enriched to GC abundances. Type II SNe expel the remain- 
ing gas and thus shut off star formation. Over several 10 8 yr, 
mass lost from massive AGB stars is returned to the gaseous 
reservoir of the GC. In addition, over this same time pe- 
riod, mass is accreted onto the GC from the ambient ISM. 
This un-processed gas is incompletely mixed with the AGB 
ejecta. All of the gas within the GC must be relatively cold 
(T ~ 10 2 K), which implies efficient cooling. However, this 
gas does not cool catastrophically for several 10 8 yr because 
of the high Lyman-Werner photon density, which photodis- 
sociates molecular hydrogen (H2). After several 10 8 yr, the 
Lyman-Werner photon density drops by more than three or- 
ders of magnitude (see Figure Q}, H2 rapidly forms, and star 
formation is triggered within the GCs gas reservoir. A sec- 
ond generation of stars is born. Type II SNe from the sec- 
ond generation clear out the remaining gao On a similar 
timescale as the formation of the second generation (several 

3 Of course, star formation must occur on a short enough timescale so 
that type II SNe do not prohibit the formation of the second generation. Star 
formation must also occur over a short timescale during the formation of the 
first generation. This has been a long-standing theoretical problem in the 
formation of bound GCs. 



10 8 yr), prompt type la SNe from the first generation begin 
to explode. These SNe then act to maintain a gas-free GC 
environment, thereby permanently ending star formation. 

Notice that GCs do not form in the centers of dark mat- 
ter halos in our scenario. Rather, in our scenario, GCs form 
wherever ISM conditions are favorable, including within gas- 
rich dwarf galaxies and interacting systems; the former being 
the presumed building blocks of the present day MW. Evi- 
dence for efficient GC formation within gas-rich dwarfs can 
be found in the high specific frequencies of GCs within dwarfs 
compared to L* galaxies (Lotz et al. 2004). 

This model contains several novel ingredients with respect 
to previous work. First, we invoke significant accretion from 
the ambien t ISM during the development of the GC gas supply 
(see a lso iPflamm-AJ tenburg & Kroupal 120091 ; iD'Ercole et alj 
2010). Second, we appeal to the importance of the high 
Lyman-Werner photon flux density in delaying star formation 
for several 10 8 yr. Third, we assert (and quantify below) that 
GC formation in the centers of dark matter halos is not a nec- 
essary condition for the formation of multiple stellar popu- 
lations because ram pressure is not effective at stripping the 
GCs gas reservoir at early times. 

We emphasize that certain aspects of this model, described 
in detail below, are not able to naturally explain the properties 
of the most massive MW GCs (e.g., wCen, NGC 2808, M22, 
and M54), which show evidence for multiple distinct stellar 
populations, and appear to require very high He abundances 
(Y ~ 0.4). We consider these objects separately in |5] 

We now provide observational and theoretical motivations 
for each of these elements. 

• As summarized in fJTJ there is no internal spread in ei- 
ther the Fe-peak elements or the a-elements (that are 
unambiguously associated with Type II SNe, e.g., Ca), 
except for the most massive GCs. This implies that the 
second generation formed out of material that was not 
enriched by the type II SNe associated with the first 
generation, nor by Type la. We conclude from this 
that type II SNe efficiently evacuated the remaining gas 
from the GC after the formation of the first genera- 
tion. Energetic arguments support this conclu sion (e.g., 
iDopita & Srnithlll986t iBaumgardt et al.ll2008t) . Notice 
that any ambient ISM that is later accreted onto the GC 
must also be unpolluted by these type II SNe. It might 
be difficult to avoid this if GCs form at the center of 
their own dark matter halos. 

• Consideration of the effects of type II SNe leads us to 
disfavor winds from massive rotating stars as the source 
of the processed material (i.e., the scenario proposed by 
Decressin et al. 2007). These stars evolve on timescales 
comparable to the onset of type II SNe. It is therefore 
very difficult to see how a second generation of stars 
could form out of the mass lost via winds from these 
massive stars both because the SNe would tend to expel 
the gas, and also because if the gas remained, it would 
be en rich ed in q-elements, which is not observed (see 
also g2]and Renzinill2008l) . We therefore favor AGB 
ejecta as the source of the processed material. 

• Evidence for accretion of pristine material comes 
from the observed abundance variations and abun - 
dance correlations within GCs. Prantz os"et alj d2007l). 
Ventura & D'Antona (2008a) and Ventura & D'Antona 
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FIG. 1 . — Evolution of the Lyman-Werner photon production rate for co- 
eval stellar populations with metallicities of log(Z/Zo ) = 0.0, — 1.0 and— 2.0. 
Rates are based on the stellar population synthesis models of Conrov et al. 
1 2009) and are normalized such that one solar mass of stars formed instantly 
at t = 0.0. Notice the rapid decline in the rate at several 10 8 yr. Lyman- 
Werner photons are capable of destroying Hi and so the precipitous drop in 
the photon production rate at t > 10 s yr implies favorable conditions for H2 
formation. 



(2009) have argued that the Na-0 anti-correlation re- 
quires dilution from prist ine material. In addition, 
iPasquini etail J2005I) and iPrantzos etafl (12007b have 
argued that the discovery of high abundances of Li and 
F in the second stellar generation requires significant 
amounts of pristine material. Their argument is based 
on the fact that these elements are fragile; Li for ex- 
ample burns at ~ 2 x 10 6 K. If the second generation 
formed purely from AGB ejecta, one might expect these 
stars to be free of Li and F. Ventura & D'Antona (2010) 
point out that Li can in fact be made in the convec- 
tive envelopes of massive AGB stars via the Cameron- 
Fowler mechanism. The implied yields are however 
strongly dependent on metallicity, stellar mass, and the 
adopted mass-loss rates. 

The need for dilution from pristine material makes it 
easier to form a second generation of stars with total 
mass comparable to the first generation, since second 
generation stars apparently do not form exclusively out 
of AGB ejecta. Notice that this un-processed mate- 
rial cannot have significantly different Fe-peak nor a- 
element abundances, which has a number of implica- 
tions. One implication is that the type II SNe that pre- 
sumably cleaned out the GC after the first generation 
was born cannot have mixed with this un-processed ma- 
terial, at least for the MW GCs, for which constraining 
data are available. 

■ Stars within GCs with abundance variations display a 
range of Na and O abundances, rather than simply 
two values as might be expected if the second gen- 
eration formed from a chemically uniform gas reser- 



voir. We interpret this result as evidence of incomplete 
mixing in the gas between the processe d and accreted 
ambient ISM material. For example, IPrantzos et al.l 
(2007) demonstrate that the full range of F and Li 
abundances can be reproduced in NGC 6742 and M4 
with a 'dilution' fac tor of ambient ISM material rang- 
ing from 0.1 to 250 IVentura & D'Antona! d2008al) and 
IVentura & D'Antona! (120091) have shown that the full 
range of the Na-O anti-correlation can be explained if 
stars form from a mixture of AGB ejecta and pristine 
material. Clearly, the GC gas reservoir must be quite 
heterogeneous prior to the second generation of star for- 
mation. 



• Over the past few years, there has been a grow- 
ing recognition that most type la SNe are prompt 
rather than delayed explosions with timescales of Gyrs 
(IScannapieco & Bildstenl 120051) . The prompt la SNe 
rate has been estimated to be roughly 10~ 12 SNe 





Mq yr ', which translates into « 10 SNe per 10 5 M, 
per 10 8 yr (e.g.. iBrandt eTaIll20Tot iMaoz et al.ll20loh . 
or ~ 10 52 erg of energy per 10 5 M© per 10 8 yr. The 
binding energy of 10 5 M© of gas with radius 1 pc — an 
extreme case — is w 10 51 erg. Once they begin to occur 
(i.e., after the delay time), prompt SNe la will therefore 
be sufficient to keep the GC gas-free. This is an im- 
portant element. If there were no prompt la SNe, it is 
not clear why a third, fourth, etc. generation of stars 
would not form from the AGB ejecta of earlier genera- 
tions. The imp ortance of prompt la S Ne has also been 
tand2008h . 



emphasized bv lD'Ercole et ; 



• Young GCs produce high flux densities of Lyman- 
Werner photons (912 < A < 1 100A). These photons are 
absorbed by H2 and ?»16% will lead to dissociation of 
the molecule. In Figure Q] we show the flux density of 
Lyman-Werner photons as a function of time for a co- 
eval stellar population. Notice the precipitous drop in 
the Lyman-Werner flux density at several 10 8 yr. We 
are interested in knowing if the flux density at < 10 8 
yr is sufficiently high to prohibit the formation of H2 in 
the bulk of the gas within GCs. This physical situation 
is analogous to photodissociation regions (PDRs) seen 
within the MW. In our case we have many, possibly 
overlapping PDRs because stars are randomly spread 
throughout the gas within young GCs. The follow- 
ing calculation is an order-of-magnitude estimate that 
proves to be sufficient for our purposes. 

We can compute the radius at which photodissociation 
of H2 is equal to the formation of H2 in a manner 
analogous to a Stromgren sphere calculation. At the 
temperatures and metallicities of interest (T w 100K, 
log(Z/Z Q ) > -3), H2 forms most efficiently on dust 
grains, with a rate coefficient of 3 x 10~ 17 ^- cm 3 s" 1 
(!Juralll975b IWolfire et al.ll2008h . The balance between 



4 The Hi formation rate and its dependence on physical parameters is quite 
uncertain. We adopt a simple linear scaling with metallicity because of the 
expected linear scaling between metallicity and dust mass fraction. This scal- 
ing is supported by Tumlinson et al. 12002) who find that the H2 formation 
rate in the LMC and SMC is a factor of ~ 10 lower than in the MW. 
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formation and destruction occurs at a radius: 
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(1) 
where Q is the photon production rate provided by 

starlight over the interval 912 < A < 1100A, n is the 
gas density, and Z is the metallicity of the gas. The 
value of Q = 10 46 s _1 is appropriate for a stellar age of 
rs 10 7 yr (see Figure[T]i. Notice that we have effectively 
assigned each GC star an average spectrum appropriate 
for a coeval population. 

This radius is to be compared to the mean interstellar 
spacing: 



'— - <%)"&)• 



(2) 



where N* is the number of stars within a radius R'. The 
gas within a young GC will form a considerable amount 
of H2 only when R LW < I. Based on the evolution in 
the flux density shown in Figure [T] this will occur only 
after several 10 8 yr, when Q decreases by a factor of > 
10 3 . Notice also that even if 90% of the Lyman- Werner 
photons are absorbed by dust (so that the effective Q is 
reduced by a factor of 10), the photon flux will still be 
sufficient to prevent significant formation of H2 during 
the early evolution of GCs. The drop in the Lyman- 
Werner flux density after several 10 s yr is so dramatic 
that our conclusion regarding H2 photodissociation is 
not sensitive to details. 

• The ~ 10 8 yr separation between first and second gener- 
ations is the natural timescale for four completely sepa- 
rate but potentially relevant processes: the orbital time 
within galactic environments, the main sequence life- 
time of massive (rs 4-8M©) AGB stars, the timescale 
for the onse t of prompt SNe la (e.g jBrandt et al.ll20 1 0t 
Maoz et all 1201 Oh . and the time when Lyman-Werner 
photon production drops precipitously for a coeval stel- 
lar population. The temporal separation between first 
and second generations cannot be significantly longer 
than 10 8 yr because otherwise prompt SNe la would 
clear out the GC gas, which would prevent the forma- 
tion of the second generation. The lack of an internal 
spread in Fe-peak elements provides additional support 
to this notion. The orbital timescale may be relevant be- 
cause tidal perturbations could be a trigger that sets off 
the second generation of star formation. Once Lyman- 
Werner photon production drops, H2 can form, and the 
gas will catastrophically cool and, presumably, form 
stars. 

In this section we have presented basic arguments favor- 
ing our adopted scenario. In the following section we expand 
upon the topic of the development of a GCs gas supply. 

3.2. GC gas reservoir growth and retention 

We will now explore the relative importance of ram pres- 
sure, Bondi accretion, and the sweeping up of the ambient 
ISM via simple geometric cross section, as a function of GC 
mass, M, relative velocity, V, between the ambient ISM and 
the GC, and density of the ambient medium, n, through with 



the GC is moving. In the following discussion we will con- 
sider Bondi accretion onto the cluster as a whole, not the ac- 
cretion of material onto individual stars. Our primary objec- 
tives in this section are twofold: to determine to what extent 
ram pressure prohibits the formation of a second stellar gen- 
eration, and to determine whether accretion from the ambient 
ISM can provide a significant source of gaseous material for 
a second generation of star formation. 

An object with mass M embedded withi n an ambient 
medium with density n will accrete matter at the lBondil(ll952l) 
rate: 



nig ~ 10 



/ M 



\ 10 5 M G 



cm 



V 



lOkms" 



M, 



©yr 



(3) 

where V can be interpreted as the sum in quadrature of the 
sound speed of the ambient medium and the velocity of the 
object through the medium. The accretion rate depends some- 
what on the equation of state of the ambient medium. 

Over 10 8 yr, Bondi accretion will result in a total mass ac- 
cretion of Mb- After 10 8 yr we therefore have: 
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(4) 



which, for a given M B and M results in a simple relation be- 
tween n and V. In this equation M represents the total mass 
within the GC. Below, we will approximate M by the "initial" 
GC mass (i.e., the mass before any gas is accreted). This ap- 
proximation yields a lower limit to the total mass accumulated 
via Bondi accretion. 

In addition to Bondi accretion, an object moving with ve- 
locity V through an ambient medium with density n and cross 
section A will sweep up material at a rate given by ni$ = pVA. 
The scenario we imagine here is one where a seed gas reser- 
voir is already present within the GC. This seed will then, by 
collisional forces, be able to sweep up additional gaseous ma- 
terial. This process will occur so long as the ram pressure is 
less than the restoring force provided by the potential well of 
the GC (see below). This ambient ISM sweeping results in a 
mass accretion rate of: 



m s « 1.5 x 10" 



V 



\cm 3 y V 10 2 tans 
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(5) 

where R is the characteristic radius of the gas reservoir within 
the GC. Notice the relatively strong dependence on radius. 
After 10 8 yr, a total mass of Ms will be swept-up, and the 
following relation will hold: 



Ms 



1.5 X l&Mr. 



cm 



V 



10 2 kms l J \3 pc 



R 



, (6) 



which, for a given Ms and R results in a simple relation be- 
tween n and V. 

A dense gaseous blob that is put in orbit within an ambi- 
ent gaseous medium will experience ram pressure from the 
ambient medium. If this ram pressure is stronger than the 
gravitational restoring force provided by the potential well in 
which the dense gaseous blob resides, then the dense gas will 
be stripped away. The critical condition, when ram pressure 
equals the gravitational restoring force, can be determined via 



pV 2 = 



GM r M g 
4ttR 4 



(7) 
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FIG. 2. — Comparison of the importance of ram pressure, Bondi accretion, and ambient ISM sweeping to the development of the gas reservoir within GCs 
(see Equations [4] [6] and[8). In this model, GCs begin with a gas reservoir with mass equal to 10% of the total mass (i.e., / = 0.1), have a half-mass radius of 
3 pc, and move with velocity V through an ambient ISM with average density n. Left panel shows results for a GC with mass 10 Mq; right panel is for a 
GC of mass 10 4 5 Mq . Shaded zones are regions of parameter space where ram pressure is strong enough to remove the gas from a GC. Lines show the amount 
of material accreted by the GC after 10 s yr by either Bondi accretion (dashed lines) or ambient ISM sweeping (dot-dashed lines). We also show in this figure 
the relation between V and n for the systems within which the GCs may be orbiting (dotted lines). This last relation is constructed by assuming gas-dominated 
self-gravitating systems (i.e., assuming V 2 = GM g /R) with gas masses indicated in the figure. 



where M» is the total stellar mass of the GC, M g and R are the 
mass and radius of gas within the GC, V is the orbital velocity 
of the GC, and p is the density of the ambient medium through 
which the GC moves. We will assume that M g = fM* and 
take R to be the half-mass radius of the cluster. This yields 
the following expression for V: 



V = 580 



-1/2 
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V 10 5 M G 
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(8) 

Notice the relatively strong dependence on radius. Observa- 
tionally, the half-light radii and masses of MW GCs are only 
weakly correlated. It is possible that any primordial mass- 
radius relation has been erased by subsequent dynamical evo- 
lution. Nonetheless, in the absence of any convincing evi- 
dence to the contrary, we assume herein that initial GC masses 
and radii are uncorrected. Unless stated otherwise, we will 
assume GC half-mass radii of 3 pc, which is the average half- 
mass radius o f GCs in the MW as determined from the MW 
GC catalog of Harris] (119961) . GCs in other galaxies also have 
avera ge half-light radii of 3 pc, independent of GC luminosity 



ra gl 



(e.g.. lMastersetal.ll2010l) . 

We emphasize that the relations derived above are only ap- 
proximate. Variations by factors of several can occur depend- 
ing on the detailed orbital properties, density profile of the GC 
and ambient ISM, and equation of state of the ambient ISM. 
These relations, and the discussion to follow, should be inter- 
preted as a guide to the trends and order of magnitude effects 
to be expected. 



The relative importance of these three processes is explored 
in Figures [2] [3] and [4] In Figure [2] we compare the impor- 
tance of ram pressure, Bondi accretion, and ISM sweeping as 
a function of the velocity, V, and ambient ISM density, n. Re- 
call that for Bondi accretion V can be interpreted as the sum in 
quadrature of the GC gas sound speed and the relative motion 
of the GC through the ambient ISM; for other mechanisms V 
should be interpreted simply as the relative velocity between 
the GC and the ambient medium. For the ram pressure calcu- 
lations, we have assumed that the mass of the GC gas reservoir 
is 10% of the total GC mass (i.e., / = 0.1; this assumption is 
discussed at the end of the section). We have also assumed a 
GC radius of 3 pc, as per the discussion above. Results in Fig- 
ure|2]are shown for GC masses of 10 4 5 M Q and 10 5 5 M Q . The 
lower mass corresponds to the lowest MW GC mass in which 
multiple populations have been found (Carretta et al. 2010c). 
Also, in Figure |2] we show the relation between n and V for 
a gas-dominated, self-gravitating system, for two gas masses 
labeled in the figure. This relation is included as a rough guide 
to the properties of the host systems in which GCs may have 
been found at their formation epoch. 

In Figure [3] we compare the same processes now as a 
function of V and GC mass, for a fixed number density of 
n = 1 cm" 3 and n= 100 cm" 3 . Finally, in Figure|4]we show the 
ratio of total accreted mass to GC mass as a function of GC 
mass at four different environments. 

There are several important conclusions to be drawn from 
these figures. First, for the present MW circular velocity 
(V w 220 km s" 1 ) ram pressure will be sufficient to remove the 
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FIG. 3. — Comparison of the importance of ram pressure, Bondi accretion, and ambient ISM sweeping to the development of the gas reservoir within GCs. This 
figure is similar to Figure|2]except for a fixed ambient density of n = 1 cm~ 3 {left panel) and n = 100 cirT 3 (right panel), and a range of GC masses, Mqc- Lines 
and shaded regions are as in Figure [2] 



gas within the GC when the GC crosses the MW disk (where 
n ~ 1 cm" 3 ) for GC masses < 1O 45 M . The circul ar velocity 
of the LMC is a factor of w 3 lower than the MW dKim et all 
1998), and so ram pressure will be able to remove the GC gas 
supply only in clusters a factor of w 3 less massivq^, or for 
M < 1O 4 M0. We therefore expect that young clusters form- 
ing with masses < 10 4 5 M in the MW should be truly co- 
eval, while young clusters of comparable mass forming in the 
LMC may host multiple populations. The critical mass for 
MW GCs is considerably higher than typically assumed be- 
cause of the high GC gas mass fraction assumed herein (10%), 
which is appropriate when considering the early development 
of GCs, when mass-loss rates are high. When one consid- 
ers ram pressure effects in present day GCs, the gas mass 
fractions considered are typically less than 1%, and s o ram 
pressure is much more effective (e.g., Frank & Gisler 1976; 
iGnedin et alj|2002h . 

Of course, the ancient MW GCs almost certainly formed in 
environments very different from the present day MW. While 
the present MW velocity is V «220kms _I , the velocity of the 
typical progenitor system in which the MW GCs formed was 
likely at least a factor of 10 lower. Thus, less massive ancient 
GCs may have also survived the effects of ram pressure (if we 



assume V, 



proj 

i3.5a 



20 kms and n = 1 cm the critical mass be- 
comes 10 j:> Mq). It is because of this fact that we emphasized 
that the critical mass of ~ 10 4,5 M Q should be manifest in clus- 
ters forming in the MW and LMC at the present epoch. With- 
out knowing in detail the properties of the progenitor systems 
in which MW GCs formed, one cannot discuss quantitatively 

5 Notice that ram pressure scales as V 2 and the restoring force scales as 
M*M g oc M \ . The critical GC mass therefore depends linearly on the veloc- 
ity. 



the effects of ram pressure on the early development of the 
ancient MW GCs. 

As mentioned in <JT] the only class of star clusters known 
not to contain multiple stellar populations are the open clus - 
ters in the MW tide Silva et al.ll2009l: iMartell & Smithll2009h . 
This population is much less massive than either the MW GCs 
or the LMC clusters, with typical masses of ~ 10 3 M Q . Based 
on FigurefJ] it is not surprising that such low mass systems do 
not contain multiple populations because ram pressure would 
be effective at stripping gas within the GC for any plausible 
formation environment. 

The right panel of Figure [3] shows the effects of ram pres- 
sure, Bondi accretion, and ISM sweeping when the ambi- 
ent ISM density is n = 100 cm -3 . These densities are very 
high compared to typical ISM environments in the local Uni- 
verse, but may be common at the epoch of massive GC for- 
mation. For example, the bulge of the MW likely formed 
rapidly at z > 2 (IZoccali et alJl2006t iBallero et alJl2007l) . has 
a stellar mass of 10 10 M^ and a half-mass radius of « 1 kpc 
tiBinnev et alj|1997h . If all the bulge stars were converted into 
gas, the gas density would be n ~ 100 cm" 3 at the epoch of 
GC formation. 

Finally, Figure [4] shows the total mass accreted within 10 8 
yrs divided by the initial GC mass as a function of GC mass. 
Depending on the formation environment, clear trends are 
predicted for the total accreted mass as a function of initial 
GC mass. In particular, if the relative velocities were low 
and the ambient densities high (upper left panel) at the epoch 
of formation, one would expect an increasing fraction of ac- 
creted material as GC mass increases. 

In this section we have assumed that the gas mass within the 
GC is only / = 10% of the GCs stellar mass. Of course, by 
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FIG. 4. — Comparison of the importance of ram pressure, Bondi accretion, 
and ambient ISM sweeping to the development of the gas reservoir within 
GCs. This figure is similar to Figures [2] and [5] except that here we show 
the fraction of mass accreted after 10 s yr as a function of GC mass for a 
variety of environments. Shaded zones are regions of parameter space where 
ram pressure is sufficient to remove the GC gas supply. Solid lines show the 
combined effects of Bondi accretion and ISM sweeping. 



the time the GC forms its second generation of stars, the gas 
mass fraction must be of order unity (or significantly greater 
than unity, depending on the star formation efficiency) in or- 
der to satisfy the observational constraint that the first and sec- 
ond generations are approximately equal in mass. Increasing 
/ to unity results in a stronger restoring force and therefore 
a reduced ability of ram pressure to remove gas within the 
GC. The critical zone for ram pressure stripping is therefore 
time-dependent. A more thorough investigation that what has 
been presented herein will therefore require time-dependent 
model calculations. In particular, the initial development of a 
gaseous reservoir will depend on the detailed hydrodynamical 
interaction between the ambient medium and AGB outflows 
in the presence of the GC potential well. 

4. CONSTRAINTS FROM INTERMEDIATE-AGE CLUSTERS IN THE 

LMC 

In this section we present new evidence for an LMC cluster 
mass threshold below which LMC clusters appear to be truly 
coeval. This result is then discussed in the context of early 
GC evolutionary scenarios. 

Multiple stellar populations have been detected in 
intermediate-a ge (~ 1 Gyr) and old LMC clusters 
(Mackev et al. 2008; Goudfr ooiTeFall 120091; iMilone et all 
2009; Mucciarelli et al. 2009). The old clusters have masses 
of « 10 5 - 5 Af Q dMackey & Gilmore 2003)); the masses of the 
intermediate-age clusters will be derived below. Intermediate- 
age clusters possess the property that age differences of ~ 10 8 
yr are readily observable in the main sequence turn-off point 
because at these ages the turn-off point is a strong function of 
time. 
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FIG. 5. — Internal age spread in intermediate-age LMC clusters a s a func- 
tion of cluster stellar mass. The age spreads are adopted from Mil one et alj 
12009) and are based on isochrone fits to ffSr-based CMDs. Stellar masses 
are derived herein. Upper limits indicate that the cluster is consistent with a 
single-age population. Notice that below log(M» /Mq )« 4.0 all LMC clus- 
ters are consistent with being coeval. The most massive cluster in this sample, 
NGC 1978, is also consistent with being coeval, although this cluster has a 
variety of peculiar properties; see the text for details. The lowest mass cluster 
has been shifted slightly in mass for clarity. 



The identification of multiple populations in intermediate- 
age clusters is important for several reasons. First, the birth 
environment of such clusters cannot be dramatically differ- 
ent from the present day LMC, in contrast with the old clus- 
ters. This opens the possibility of making more direct links 
between the clusters and their environment than is possible 
with the old clusters. Second, as discussed in ^2.1.11 these ob- 
servations argue strongly against GC formation at the center 
of their own dark matter halos as a means to produce multiple 
p opulations. 

IMilone et alj (2009) has recently analyzed archival HST im- 
ages of 16 intermediate-age clusters in the LMC. These au- 
thors derive not only average ages, metallicities, distances, 
and reddening values, but also the internal spread in age for 
each cluster in their sample. Only five of the clusters in their 
sample were consistent with being coeval. We are interested 
in knowing whether the internal age spreads measured by 
Milone et al. correlate with the mass of the cluster. We have 
therefore derived stellar masses for the clusters in Milone et 
al., in the following way. The average age of each cluster is 
converted into a V -band mass-to-lig ht ratio, M/Ly, via the 
stellar populatio n synthesis model s of IConrov et al.1 (120091) as 
updated in Conrov & Gun n (2010). We assume Z = 0.006 for 
the young clusters in the LMC and a Kroupa (2001) IMF. We 
then adopt V-band magnitud es for the Milone et a l. clusters 
from the data compilation of Ivan den Berghl dl981l) (there are 
fourteen clusters in common). These photometry are then cor- 
rected for extinction and converted to absolute luminosities 
via the E(B - V) values and distance moduli listed in Milone 
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et al. Finally, stellar masses are estimated by combining the 
expected M/Ly and absolute magnitudes. The majority of 
the clusters have similar ages and so the estimated M/Lys do 
not vary considerably across the sample. The following result 
therefore does not depend critically on the adopted modeling 
approach. 

In Figure we show the relation between th e total stellar 
mass and internal age spread for the clusters in Milo ne et al.l 
(2009). It is remarkable that all clusters with masses < 10 4 M Q 
show no evidence for multiple populations. In the LMC, there 
appears to be a critical mass below which the formation of 
multiple populations is suppressed. This critical mass is fully 
consistent with our estimates of the effects of ram pressure 
effects in the LMC, where for n = 1 cm" 3 and V = 70 kms" 1 
we find a critical mass of M sa 10 4 M Q . Given the sensitivity 
to parameters such as the GC radius, a more quantitative com- 
parison must await detailed observations of these clusters. 

The most massive cluster in Figure [5] NGC 1978, also 
sho ws no evidence for m ultiple populations. It is worth noting 
that lMilone et al.1 (120091) extract the CMD for their target clus- 
ters from the cluster core for all clusters except for NGC 1978 
because of crowding. Since the second generation appears to 
be more spatially concentrated than the first (see EjTJ, this may 
explain the observed lack of significant internal age spread in 
this cluster. This cluster ha s a num ber of peculiar properties, 
including a high ellipticity (Geisler & Hodge 1980), a lack of 
S stars, and odd carbon and oxygen ratios in its AGB stars 
dLederer et al.ll2009t) . Owing to all of these facts, we are re- 
luctant to place much weight on the result that this massive 
cl uster does not show obviou s signs of internal age spreads. 

iBastian & de Minkl (120091) raise the possibility that the ob- 
served broadening of the CMD of intermediate-age LMC 
clusters is in fact due to stellar rotation, and not multi- 
ple populations. We believe this scenario to be implausible 
for several reasons. First, our results in Figure [5] strongly 
suggest that the width of the main sequence turn-off is re- 
lated to the total cluster mass. This correlation is difficult 
to explain if the width was due to stellar rotation. Sec- 
ond, no such broadening is observed in the intermediate-age 
MW open clusters, including the Hyades, Pleiades, and Prae- 
sepe (e.g., Vandenberg & Bridges 1984; G riffin et al.l 1 19881: 
Pinsonneault et al. 2004; An et al. 2007). These MW open 
clusters have masses of order 10 3 M Q , and so in our scenario 
they should be truly coeval, as observed. The scenario of Bas- 
tian & de Mink can be reconciled with the open cluster data 
and low mass LMC data only by supposing an ad hoc cor- 
relation between stellar rotation rate and cluster mass, which 
seems implausible. 

5. SUMMARY, OPEN ISSUES AND FUTURE DIRECTIONS 

In this work we have presented a comprehensive model for 
the early (several 10 s yr) evolution of massive star clusters. 
Our model considers the importance of type II and prompt 
type la SNe, accretion onto the GC from the ambient ISM, the 
ability of a GC to retain its internal gas supply in the face of 
ram pressure, and the effect of the Lyman- Werner photon flux 
density on the ability of the young GC gas to form molecules 
and, ultimately, stars. 

This model definitively addresses two of the three major is- 
sues in early GC evolution identified in £|2] how GCs can re- 
tain a gaseous reservoir in the face of ram pressure, and which 
stars are responsible for processing material to T > 10 7 K. 
From consideration of the formation environments of both old 



MW GCs and intermediate-age LMC clusters, we have shown 
that ram pressure stripping naturally explains the observed bi- 
furcation between clusters that do and do not show evidence 
for multiple stellar populations. Based on several independent 
arguments, including the similar ~ 10 8 yr timescale of many 
physical mechanisms and the lack of internal spread in Fe and 
Ca in most clusters, we strongly favor massive AGB stars as 
the source of the processed material. 

The final open issue is in understanding how a second gen- 
eration can form with a current total mass comparable to the 
first generation, at least for old MW GCs where data are avail- 
able. We have considered accretion from the ambient ISM as 
a viable mechanism to provide copious amounts of gas to the 
GC. However, our solution to this last issue currently has lit- 
tle direct empirical verification, and may in fact have trouble 
reproducing the observed correlations in abundance ratios, as 
discussed below. We now discuss a variety of testable im- 
plications of our proposed solution, and comment briefly on 
several open issues. 

As a consequence of significant accretion from the ambi- 
ent ISM, we expect several trends with GC mass (see e.g., 
Figure |4). For example, as the fraction of pristine gaseous 
material increases, we expect a shorter Na-O correlation, a 
greater abundance of Li and F, and in general we expect the 
stars with anomalous abundances to be less anomalous. The 
sign of the trend with mass depends on the dominant accre- 
tion process, which in turn depends on the formation environ- 
ment. Bondi accretion scales as M 2 while the mass provided 
by AGB winds scales with M so the pristine fraction increases 
with mass. ISM sweeping scales as R 2 which is only weakly 
correlated with M at the present epoch, and so if this process 
is dominant we would expect the pristine fraction to decrease 
with mass. 

Observations of young and intermediate-age clusters will 
allow identification of the relevant mechanism because the 
formation environment of young clusters is not so dissim- 
ilar from their present environment. The old GCs within 
dwarf spheroidals and dwarf irregulars may also shed light 
on this issue. The GCs within dwarf galaxies almost certainly 
formed where they are now observed since the stellar accre- 
tion/merger rate onto dwarfs is expected to be very low in a 
ACDM cosmology, and so the present day conditions of the 
host dwarf galaxies cannot be very different from the forma- 
tion environment of the GCs. 

Carretta (2006) and ICarretta etaP (1201 Obi) have investi- 
gated the extent of the Na-O correlation as a function of global 
parameters and have found evidence that the correlation is 
more extended in higher mass GCs, and for GCs with more 
extended orbits. This result is consistent with ISM sweeping 
being the dominant accretion process, but we caution that the 
data show large scatter, and the present GC mass is poorly 
correlated with its mass at formation due to orbit-dependent 
mass-loss effects. In general, we expect that the extent of 
the Na-O correlation should depend on the formation envi- 
ronment at fixed GC mass, since the formation environment 
determines the amount of accreted material. Another effect 
might be the mass-dependent ability of a GC to retain the 
winds from AGB stars. In any event, it is clear that such 
trends hold the promise of isolating which physical process 
dominates the accretion of pristine material onto GCs. 

If accretion from the ambient ISM was important, then it is 
somewhat puzzling why the Fe and Ca abundances are so uni- 
form between the first and second generations. The only plau- 
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sible explanation is that at early times, during the formation 
of the ancient MW GCs, the spread in Fe and Ca abundances 
within the MW progenitor system was quite small. A clear 
prediction of the accretion scenario is that the intermediate- 
age LMC clusters showing multiple populations should show 
a much larger internal spread in Fe, Ca, etc. abundances, on 
the order of the spread in abundances of these elements within 
the LMC as a whole. 

The MW bulge GC Terza n 5 may provide additional in- 
sight. iFerraro et al.l (120091) has convincingly shown that 
this GC has a split horizontal branch (HB) with the more 
luminous branch having a much higher Fe abundance 
([Fe/Hl~ +0.3) t han th e less luminous branch ([Fe/H]~ -0.2). 
iD'Antona et alj d2010l) argue that the split HB is consistent 
with an internal age spread of several 10 8 yr and enhanced He 
in the metal-enriched branch. Since this cluster is in the MW 
bulge and Fe-rich, it probably formed in the bulge. D' Antona 
et al. suggest that the second stellar generation may have ac- 
quired its high Fe abundance via accretion from the ambient, 
Fe-rich ISM. Self-enrichment is unlikely because the required 
Fe mass to produce the observed high Fe abundance would 
require so many type II SNe that the cluster would easily be- 
come unbound, unless Terzan 5 was significantly more mas- 
sive in the past. Terzan 5 may therefore be an ideal cluster to 
look for further evidence for the importance of ambient ISM 
accretion in the formation of multiple stellar populations, e.g., 
by observing the Li and F abundances in this cluster. If the 
proposed scenario for Terzan 5 is generic then we might ex- 
pect to find significantly separated Fe abundances within the 
other bulge metal-rich GCs as well. 

In the present work we have not attempted to make spe- 
cific predictions for the extent of correlations amongst vari- 
ous elemental abundances. Such predictions would provide a 
powerful constraint on the model. Qualitative considerations 
suggest that substantial accretion from the ISM may yield a 
very short Na-0 anti-correlation, owing to the ISM having 
abundance ratios similar to the first generation stars. Unfortu- 
nately, uncertainties in the AGB yields greatly complicate any 
attempted comparison to observed abundances. A fruitful av- 
enue for future work will consist of a detailed comparison be- 
tween predicted and observed ab undance correlations i n light 
of the uncertain AGB yields Csee lD'Ercole et alE OlO. for an 
initial attempt in this direction). 

The most massive GCs, including wCen, NGC 2808, M22, 
and M54, deserve special mention. Each of these clusters 
show multiple, distinct sequences in the CMD, and the former 
two are unique in that they display multiple main sequences. 
The CMD morphology of these clusters suggests very high 
He abundances in the second (and third) stellar generations. 
The discrete sequences and high He abundances argue against 
significant dilution from ambient ISM accretion. This can be 
accommodated in our model if ISM sweeping is the domi- 
nant accretion mechanism because this mechanism becomes 
increasingly less important at high GC masses (see the bot- 
tom panels of Figure HJ. Although we are then left with the 
original problem of explaining how so many second genera- 
tion stars could form from the mass lost by the first genera- 
tion. We speculate that perhaps the most massive GCs formed 
at the centers of their own dark matter halos, and contained 



even more stars in their past. The precursors of these systems 
could be nuclear star clusters, which are common in low-mass 
galaxies. Regardless of these details, we stress caution when 
attempting to incorporate the massive GCs into any frame- 
work for the early evolution of GCs since many of their ob- 
servational characteristics differ qualitatively from the lower 
mass GCs. 

An exciting direction for future observational work is 
characterizing the incidence of multiple populations within 
the numerous young and intermediate-age clusters in M31. 
These clusters span a range in mass from 10 3 M Q to 10 5 M Q 
(Caldwell et al. 2009). Based on ram pressure arguments, we 
expect the higher mass clusters to show evidence for multi- 
ple populations, but not the lower mass clusters. Assuming 
n = 1 cm" 3 and V = 250 km s" 1 for disk of M31 dWidrow et alj 
2003), we predict that the critical mass will be ss 1O 4 ' 5 M . 
However, as stressed in previous sections, the relevant veloc- 
ity is the relative velocity between the ambient ISM and the 
GC. Since many of the young M31 GCs are orbiting within 
the disk of M3 1 , their relative velocity will be considerably 
less than the circular velocity, and so the critical mass may be 
considerably lower. Confrontation of our model predictions 
with the properties of the M3 1 GCs must take these details 
into account. 

Since these M31 clusters have not experienced significant 
dynamical evolution, we may expect stronger correlations 
between the extent of any abundance anomalies and cluster 
mass. For the youngest clusters, with ages of several 10 8 yr, 
we might even hope to directly observe the formation of the 
second generation. Such an observation would provide defini- 
tive proof that GCs need not form within dark matter halos to 
produce multiple stellar generations, and would also demon- 
strate that AGB stars are the polluters, since the AGB pol- 
luter scenario is expected to produce a second generation on 
a timescale of several 10 8 yr. 

In the past several years it has become clear that star clus- 
ters, once thought to be simple systems, in fact show an in- 
ternal complexity that increases with increasing mass. Obser- 
vations indicate that the lowest mass systems, i.e., the open 
clusters, appear to be truly coeval and mono-metallic. At 
higher masses one observes age spreads of several 10 8 yr 
and internal variation in the light elements. At still higher 
masses (> 10 6 M©) it appears that the cluster potential well is 
deep enough to retain SNe ejecta and hence self-enrich, per- 
haps because these massive clusters form at the center of their 
own dark matter halos. This last phenomenon is observed not 
only in the most massive MW GCs, but also appears to be 
present in th e most massive GCs within other galaxies as well 
(IStrader & S mith 2008; Bailin & Harris! 2009). At the very 
least, it is now abundantly clear that star clusters are not sim- 
ple systems. 
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and thank Alvio Renzini for comments on an earlier draft. The 
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archive at arXiv . org. 
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